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1.  INTRODUCTION 

As  part  of  a  study  of  acoustic  processing  techniques  it  was  necessary  to  sake 
estimates  from  experimental  data  of  the  ratio  of  received  signal  power  to  noise 
power  per  hertz.  One  method  is  to  record  at  a  sensor  a  sinusoidal  tone  emitted 
by  a  distant  source.  The  received  signal  is  then  digitised  and  the  power 
spectrum  confuted  from  a  set  of  data  samples.  The  sample  rate  and  the  number 
of  data  samples  are  chosen  so  that  an  appropriate  frequency  bin  of  the  spectrum 
will  contain  the  total  signal  power  apart  from  side  lobes  of  negligible  height. 
This  signal  bin  will  also  contain  a  component  of  noise  power.  Two  alternatives 
are  now  possible: 

First,  an  estimate  (z)  of  the  ratio  of  signal  to  noise-per-hertz  can  be  obtained 
from  the  single  spectrum  by  evaluating 


z 


(1) 


or 


z 


(2) 


where  s2  is  the  power  including  noise  measured  in  the  signal  bin. 

u  is  the  average  power  in  frequency  bins  near  the  signal  bin  but  not 
containing  extraneous  signals. 

i  is  the  bin  width  in  hertz, 
r 

Secondly,  a  number  of  spectra  from  consecutive  sets  of  data  samples  can  be 
obtained.  An  average  spectrum  can  then  be  formed  and  an  estimate  of  z  can  be 
obtained  from  using  either  equation  (1)  or  (2). 

In  this  paper  a  probabilistic  analysis  is  presented  for  the  first  alternative. 
A  number  of  assumptions  were  found  to  be  necessary  to  arrive  at  a  satisfactory 
result; 

(1)  The  received  signal  tone  has  constant  amplitude  during  the  time  over  which 
an  estimate  of  signal-to-noise  ratio  is  made. 

(2)  The  number  of  spectral  bins  used  for  spectral  analysis  is  large, 
e.g.  the  number  of  samples  used  is  >  1024. 

(3)  The  noise  field  as  recorded  by  the  sensor  is  isotropic,  stationary, 
continuous  and  white. 


2.  TOE  DISTRIBUTION  OF  z  WHEN  ESTIMATED  FROM  A  SINGLE  SPECTRUM 

The  power  spectrum  of  a  sampled  time  series  is  the  output  from  a  digital 
algorithm  and  consists  of  a  tabular  sequence  of  power  estimates  for  each 
frequency  bin. 

Following  Rice(ref.l)  the  probability  density  function  of  the  envelope  of  a 
signal  tone  in  narrow  band  noise  is  given  by 

fg(x)  -  2Alx  exp(-A2(xJ+K,J)  I0(2AJKx),  x  >  0  (3) 
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where  K  is  the  amplitude  of  the  signal  tone 
£^2  is  the  mean  square  noise  amplitude 
1q  is  the  modified  Bessel  function  of  zero  order. 

The  noise  can  be  assumed  to  have  zero  mean  and  hence  the  noise  variance  is  - 

2i 

The  density  function  for  the  envelope  of  narrow  band  noise  is  obtained  by 
putting  K  equal  to  zero  in  equation  (3),  i.e. 


f^(x)  =  2A2x  exp(-A2x). 


The  estimator  for  signal-to-noise  ratio  is 


z 


(4) 


(5) 


where  s2  is  the  power  in  the  frequency  bin  containing  the  signal  tone 

n2  is  the  power  in  any  bin  that  contains  noise  only 

X 

p  is  the  number  of  bins  used  to  estimate  the  average  noise  power 

—  is  the  bin  width  in  hertz, 
r 

Therefore,  the  estimator  (z)  is  the  ratio  of  the  signal  tone  power  to  noise- 
power-in-a-1  Hz-band. 

As  the  bins  have  a  width  that  is  generally  small  in  comparison  with  the  signal 

bin  centre  frequency,  equations  (3)  and  (4)  can  be  used  to  derive  the  distributions 

of  s2  and  nj^  in  equation  (5) . 

Generally,  for  probability  density  functions  where  fx(x)  =  0  for  x  <  0 

fX2W  =  2^fXCx^>  W 

This  can  be  found  in  Parzen(ref . 3) . 

Hence 

fgj  (x)  =  A2  exp(-A2  (x  +  K2 ))  I^A2^)  (7) 


and 


*Ni(x)  =  A2  exp(-A2x) .  (8) 

i 


Equation  (8)  is  the  density  of  a  function  distributed  as 
In  Appendix  I  it  is  shown  that  u,  where 


P 


i*l 


JO. 


-  3  - 
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is  distributed  as  X^p  i/(A2 V*)  ’  i,e' 


fu(x)  = 


xP-1  exp 


(10) 


The  probability  density  function  for  the  ratio  s2/u  can  now  be  obtained  by 
applying  the  classical  relationship  given  in  Parzen  for  the  density  of  the  ratio 
of  2  independent  positive  random  variables 


f^yCy)  =  f  x  fx(yx)  fy  (x)  dx 
o 


(11) 


i.e. 


C 

fs*/u(y)  s  J 


dx  x  A2  exp(-A2(xy  K2))  I0(2A2KV$cy) 


A2PnP  n-1 

X  xp  1  exp(-A2px) 


(12: 


A2p+2pp  exp(-A2K2)  ,« 

- /  dx  xp  exp(-A2x(y+p)) 


Iq(2A2K  V^y)  ♦ 


(i3: 


For  ease  in  writing  put 


a  = 


A2p*2pP  exp (-A2  K2 ) 
T(P) 


(14 


a 

y 


A2  (y  +  p) 
2A2  K  y/y. 


(15 

(16 


Then 


fS2  /U  *-y  ^ 


<X 

“/ 

"  n 


dx  xp  exp(-0x)  In(7  ^5t) 


(17 


Making  the  substitution  u  =  >/x,  2u  du  =  dx 
then 


fs’/u<y)  ■  2" 


OO 

f  du  u2p+1  exp(-0u2 )  IqCTu)  • 
J  o 


(18 
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Using  (9.6.3)  of  Abramowitz  and  Stegun(ref .2) ,  this  equation  can  be  rewritten  as 


fgi/uCy)  =2 a  f  du  u2p+1  exp(-^u2)  JQ(i7u).  (19) 

J  o 

The  integral  can  be  simplified  and  evaluated  by  replacing  the  Bessel  function 
with  its  series  expansion  and  reversing  the  order  of  integration  and  summation. 
The  resulting  integral  in  the  general  term  can  be  evaluated  using  the  integration 
formula  (7.4.5)  in  Abramowitz  and  Stegun.  Hence,  after  some  rearrangement 


OO 


j=0 


Rewriting  the  Gamma  function  in  terms  of  Pochhammers  symbol  defined  for  example  in 
(6.1.22)  of  Abramowitz  and  Stegun  gives 


i.e. 


(21) 


(22) 


where  M(.,  .,  .)  is  the  Confluent  Hypergeometric  function  known  as  Rummer's 
Function  and  whose  series  expansion  is  given  for  example  in  (13.1.2)  of 
Abramowitz  and  Stegun.  Substitution  from  equations  (14),  (15)  and  (16)  now 
enables  the  density  function  to  be  written  in  terms  of  its  primary  parameters  as 


fs'/u<y)  -  “jm  h(p*i.  1.  7^)-  (25) 

The  distribution  function  for  z  now  follows  from  the  relation 

Wy>  ■  ifx(5r) 

i.e. 

tzU)  '  („  !vl  7)^  “p(‘a!k,)  m(p*‘-  *•  (24) 

where  <  z  <  °°. 
r 

Typical  shapes  taken  by  this  distribution  are  shown  in  figure  1.  These  were 
derived  using  the  simple  program  given  in  Appendix  II. 
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3.  THE  MEAN  AND  VARIANCE  OF  z 

The  mean  and  variance  of  z  are  derived  readily  from  the  mean  and  second 
moment  of  S2/U.  The  moments  of  S2/U  are 


E(yn)  =  [  dy  yn  -P-^J  exp(-A2K2)  m(P*1,  1. 

J  n  (y+pr  v  7 


(25) 


Substitution  of  the  series  exjjansion  for  M(.,  . ,  .)  and  altering  the  order  of 
integration  and  summation  gives  rise  to  an  integral  in  the  general  summation  term 
which,  on  evaluation,  forms  a  0 -function. 

The  0-function  can  then  be  written  in  terms  of  T-functions.  After  some 
algebraic  maniPulation : 

E(yn)  =  Pn+1  exP(-A2K2)  0(P-n,n+l)  M(n+1,  1,  A2K2).  (26) 

The  details  of  the  derivation  are  given  in  APPendix  III. 

The  mean  of  y  is  now 

E(y)  =  p2  ex]3(-A2K2)  0  (P- 1 ,2)  M(2,  1,  A2K2).  (27) 

From  (13.4.1)  and  (13.6.12)  in  Abramowitz  and  Stegun 

M(2,  1,  A2K2)  '=  (A2K2+1)  exP (A2 K2 ) . 


From  (6.2.2)  in  Abramowitz  and  Stegun 


0 (P~ 1,2)  = 


^(p-1)  T(2) 
T(p+1) 


1 

P(P-1)‘ 


Hence  after  some  rearranging 


E(y) 


(A2K2+1)P 
(P-1)  * 


The  second  moment,  derived  in  a  similar  fashion,  is 

E(y2 )  =  P3  exp(-A3  K2)  0(P-2,3)  M(3,1,A2K2). 

From  (13.4.1)  and  (13.6.12)  in  Abramowitz  and  Stegun 

M(3, 1,A2 K2 )  =  *s(2  +  4A2 K2  ♦  A4K*)  exP(A2K2). 


From  (6.2.2)  in  Abramowitz  and  Stegun 


P (P-2, 3) 


rq>-2)  T(3) 

T(p+1) 


2 

P (P-1) (P-2)  * 


(28) 


(29) 
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After  soae  rearranging 


E(y*) 


p2  (2  ♦  4A*K2  ♦  A4**) 

— (p-i)  (p-23 - 


(30) 


The  aean  and  second  aoaent  of  z  now  follow  froa  application  of  the  properties  of 
the  expection  operator  E  to  the  equation 


z 


i 


.e. 


E(z)  =  7  E(y)  -  i 

p{l  ♦  A2K2)  1 

r(p-l)  '  r 


froa  equation  (28) 


Again 


E(z2 )  *  p  Ely2)  -  p  E(y)  +  p  . 


The  variance  of  z  is  then 


o2  (z)  =  E(z2)  -  (E(z))2 

=  p  E (y2)  -  p  E (y)  ♦  p  -  p  (E(y))2 

♦  p  E(y)  -  ^ 


(31) 


i.e. 


*’(*)  *  p  Efy2)  -  p  (E(y))2 . 

Substitution  froa  equations  (28)  and  (30)  gives  after  sane  rearrange 


o2 (z)  =  e^+.^p.**4*4)  . 

r2  (p  -  l)2  (p  -  2) 


(32) 
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4.  TOE  LIMITING  DISTRIBUTION  OF  z  IN  TOE  ABSENCE  OF  A  SIGNAL  TONE 

When  the  signal  has  zero  amplitude,  the  signal-to-noise  ratio  is  zero  and  the 
probability  density  function  and  associated  mean  and  variance  are  simply  obtained 
by  putting  A2K2  to  zero  in  equations  (24),  (31)  and  (32)  respectively. 

Hence  the  probability  density  function  has  the  form 


fz(z) 


,P+1 


(rz  +  p  +  1) 


p+1 


(33) 


where  <  z  <  °°. 
r 

The  mean  is  given  by 


E(z) 


1 

r  (p  -  1)* 


and 


(34) 


the  variance  is  given  by 


o1  (z) 


_ p! _ 

r2  (P  -  1)2(P  -  2)' 


(35) 


5.  TOE  LIMITING  DISTRIBUTION  OF  z  WHEN  TOE  INPUT  SIGNAL-TO-NOISE 

RATIO  IS  LARGE 

For  large  signal-to-noise  ratio  the  relation  (13.1.4)  in  Abamowitz  and  Stegun 
can  be  used  to  rewrite  f^(z)  as  expressed  by  equation  (24)  to  give 


JP+1 


exp(-A2  K2 )  -L  <>gllL  expj^I^] 


(rz  +  p  +  1) 


p+1 


(p+1)  (rz+p+l)p  L  (rz+p+1) 


i 


i.e. 


fz(Z) 


pP+1(rz+l)P_  )  -A 


2p+l  exp  (_  rz+p+ 


2  K2  P  ] 

WJ 


(&>) 


pj (rz+p+1) 

The  mean  and  variance  follow  simply  from  equations  (31)  and  (32),  i.e. 

EU)  =  P^L 

r(p-l) 


(37) 


and 


o2 (z)  = 


52A4K4 


r2  (p-1)2  (p-2) 


(38) 


mmm 
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6.  DESIGN  OF  EXPERIMENTS 

In  order  to  use  the  preceding  analysis  in  the  design  of  experiments,  the 
coefficient  of  variation  of  the  distribution  (equation  (24))  should  be  examined 
for  selected  values  of  r  and  p. 

The  coefficient  of  variation  is  defined  as 


C  *  E(z) 

and  is  a  measure  of  the  dispersion  of  a  distribution  function  when  the  expected 
value  is  the  unit  of  measurement. 

From  equations  (31)  and  (32) 

c  =  P  ( P  »  2PA*  K2  +  A4  K4 

pA2  K2  +  1  V  p  -  2 

A  typical  behaviour  of  C  is  shown  in  figure  2.  It  would  appear  that  any 
improvement  in  handling  low  signal-to-noise  ratio  by  using  higher  values  of  p 
is  not  possible. 

An  alternative  way  of  looking  at  the  design  of  experiments  is  shown  in 
figure  3.  Here  the  bias  in  the  mean  of  the  measurement  is  clearly  seen.  Also, 
the  problem  in  the  long  right-hand  tail  of  the  distribution  giving  an  increased 
chance  of  obtaining  negative  values  for  the  estimate  of  signal-to-noise  ratio  is 
evident.  This  latter  problem  is  also  illustrated  by  equation  (34)  and  (35) 
where  E(z)  -  a  is  always  negative  and  equations  (37)  and  (35)  where  E(z)  -a  is 
always  positive. 


7.  CONCLUSION 

An  analysis  has  been  presented  of  a  frequently  used  estimator  for  signal-to- 
noise  ratio.  The  result  in  closed  form  is  new  as  far  as  the  author  can 

determine.  A  simple  technique  for  examining  any  proposed  experiment  is  then 
described  and  points  to  the  inability  to  measure  signal-to-noise  ratio  with  any 
precision  at  low  values  of  the  input  signal-to-noise  ratio.  This  points  to 
the  necessity  of  examining  the  use  of  multiple  estimates  of  the  power  spectrum 
to  obtain  an  average  spectrum  from  which  an  estimate  of  signal-to-noise  ratio  can 
be  made. 

Some  progress  has  already  been  made  with  this  problem  and  it  is  proposed  to 
publish  a  paper  on  this  topic  shortly. 


t 
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APPENDIX  I 


THE  DISTRIBUTION  OF 


If  x  is  distributed  as  x2  then  x  has  a  probability  density  function 

n,a 


fY(x)  = 


‘XW  '  2n72  on  r(n/2) 

Hence  0  has  a  probability  density  function 

£*  e*p{#] 


x"72-1  exp  ran. 


fY(x)  = 


The  addition  of  p  samples  from  x^  c  is  equivalent  to  adding  2p  squared  terms  where 
each  term  is  a  sample  from  NCO,  o );  i.e.  the  sum  of  p  samples  from  x^  g  is 
distributed  as  xl  n  and  the  probability  density  function  of  the  sum  is  given  by 


f*(x)  -  2P  „*>  r(p) 


1  “»[»»] 


The  probability  density  function  of  the  average  is  found  using  the  general 
relation  (Parzen) 

faX«  *  a  fX  <!> 


i.e. 


f«w 


PP-‘  xp‘‘  -»(«?•) 


2V  or  r(p) 


T r<p) 


)-l 


exp 


r_-x _ 


i.e.  U  is  distributed  as  x 


2p 
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APPENDIX  II 

EVALUATION  OF  PROBABILITY  DENSITY  FUNCTION  EQUATION  (4) 

NOTE  THAT  THE  CONFLUENT  HYPERGEOMETRIC  FUNCTION  IS 
EVALUATED  RECURSIVELY  IN  PROCEDURE  FUNC  USING 
EQUATION  (13.4.1)  IN  HANDBOOK  OF  MATHEMATICAL 
FUNCTIONS  BY  ABRAMOWITZ  AND  STEGUN,  DOVER  PRESS,  1970 

SNRPDF :  PROCEDURE  OPTIONS  (MAIN) : 

Z  =  -1; 

R  =  1;  /*  1/R  *  FREQUENCY  BIN  WIDTH 

IP  *  18;  /*  IP  =  NO  OF  BINS  OVER  WHICH  NOISE  IS  AVERAGED 

Y  =  2;  /*  Y  =  INPUT  SNR  IN  DECIBELS 

K  =  10** (Y/10) ; 

DO  I  =  1  TO  50; 

A  =  K*(R*Z+1)/(R*Z+IP+1) ; 

B  =  FUNC (IP, 1, A) ; 

C  =  ((IP/ (R*Z+IP+1))**(IP+1))*EXP(-K)*B; 

PUT  EDIT  (Z,C) (SKIP(l) ,F(2) ,X(3) ,F(12,6)) ; 

Z  =  Z  ♦  1;  END;  END; 

FUNC,  PROCEDURE  (IP,1,X); 

Ml  =  EXP(X) ;  M2  =  (2+X)*EXP(X) ; 

DO  I  =  3  TO  IP; 

M3  =  ((2*(I-1)-1+X)/ (I-1)*M2  -  ( 1- (1/ (I-1)))*M1; 

Ml  =  M2;  M2  =  M3;  END; 

RETURN (M3);  END  FUNC; 

END; 


i 


WSRL-0139-TM 


12  - 


WSRL-0139-TM 
Figure  1 


Figure  1.  Probability  density  functions  for  r  =  1,  p  =  18  and 
input  signal-to-noise  ratios  of  0  dB,  3  dB,  6  dB 
and  10  dB 
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